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Numerical and theoretical studies of a one-dimensional spin model with locally broken spin sym- 
metry are presented. The multi-species annihilating random walk transition found at zero branching 
rate previously is investigated now concerning the cluster behaviour of the underlying spins. Generic 
power law behaviors are found, besides the phase transition point, also in the active phase with ful- 
fillment of the hyperscaling law. On the other hand scaling laws connecting bulk- and cluster 
exponents are broken - a possibility in no contradiction with basic scaling assumptions because of 
the missing absorbing phase. 
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INTRODUCTION 



The study of non-equilibrium model systems has at- 
tracted great attention in recent years. A variety of phase 
transitions have been found characterized by critical ex- 
ponents, both static and dynamic. Of special interest 
are transitions from a fluctuating active state into an ab- 
sorbing one. A wide range of models with transitions into 
absorbing states was found to belong to the directed per- 
colation (DP) universality class l] . Another universality 
class of interest is the so-called parity conserving (PC) 
class [1,11,^1^. The mostly studied particle model in this 
class is branching annihilating random walk (BARW) in 
one dimension (Id) with an even number of offsprings 
{2 A ^ 0, A — > 3 A, in the simplest case). The first exam- 
ple of model systems exhibiting PC-type transition was 
given, however, in two Id cellular automata by Grass- 
berger The prototype spm-model for PC-type phase 
transitions was proposed by one of the authors ^] by in- 
troducing a class of non-equilibrium kinetic Ising models 
(NEKIM) with combined spin-flip dynamics at zero 
temperature, T = and Kawasaki spin-exchange kinet- 
ics at T = oo. 

Transitions between active and absorbing phases have 
been, however, mostly studied in particle-type models. 
The N-BARW2 model is a classical stochastic system 
of N types of particles with branching annihilating ran- 
dom walk. For N > 1 N types of particles Ai perform 
diffusion, pairwise annihilation of the same species and 
branching Ai Ai + 2Aj with rate a for i ~ j and 
with rate <j' /{N — 1) for i ^ j. According to fleld theory 
[lol| in this model the rate a flows to zero under coarse- 
graining renormalization which implies that the model 
is always active except for the annihilation fixed point 
at fj =0. It forms a universality class, the so called 
N-BARW2, different from DP and PC, with well known 
bulk critical exponents in Id. 
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In the NEKIM model a global asymmetry of the spins 
(magnetic field) is known to change the PC transition 
into the DP type. [H El ■ The question arises what 
is the effect of a local breaking of the spin symmetry 
in such a spin system. The first indication in this di- 
rection has come from a work of Majumdar et al. [l^ 
who studied the coarsening dynamics of of a Glauber- 
Ising chain with strong asymmetry in the annihilation 
rate maximally favouring '-' spins (MDG- model). These 
authors found the result that while the domains still 
coarsen as t^/"^, the '-' domains coarsen slightly faster as 
t^/^log(i). As a result at late times, the system started 
from a random initial state decays into a fully compact 
state where all spins become '-' in a slow logarithmic way 
l/log(i). ^ 

In a previous paper |l4| the authors presented an 
asymmetric spin-model (NEKIMA) by generalizing the 
NEKIM model which includes as a special case the MDG 
model. In NEKIMA there is local spin-asymmetry both 
in the annihilation rate (favouring '-' spins) and the 
diffusion- like spin-fiip rate (favoring '-f ' spins) and thus 
acting oppositely. Global scaling properties of the model 
have been investigated numerically as well as using clus- 
ter MF approximation. The N-BARW2 transition, for 
which no spin model had been known previously, was 
found at zero spin-exchange rate. In the present pa- 
per we further investigate this model at such parame- 
ter values for which in the original NEKIM model PC- 
type transition takes place. In the plane of the spin- 
asymmetry parameter and kink-branching probability we 
have found, by computer simulations as well as by clus- 
ter mean-field calculations, a phase diagram showing a 
reentrant directed-percolation line. Our main purpose, 
however, has been to investigate the cluster development 
properties a) at and in the vicinity of the N-BARW2 line 
andb). in the rest of the parameter space considered. For 
the mean population size nit) ^ , for the mean square 
spreading of spins R'^{t) ^ V' and for the survival prob- 
ability P{t) ^ generic scaling behaviour has been 
found via computer simulations in (almost) the whole 
plane of the phase diagram with fulfillment of the hy- 
perscaling law. Upon crossing the line of zero branching 
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rate (where the phase transition takes place), however, 
dynamic scaling is found to be violated concerning laws 
connecting bulk exponents and cluster ones. We trace 
back such a possibility to the circumstance that the ab- 
sorbing phase is missing by the N-BARW2 transition. 

II. THE MODEL AND PREVIOUS RESULTS 

The general form of the Glauber spin-flip transition 
rate in one-dimension for spin Si sitting at site i is 
(s, =±1): 

T ~ 1 

w{si,Si-l,Si+l) = -(H-(5si_iSi+i)[l--Sj(s,;_i+Si+i)] 

(1) 

at zero temperature. (Usually the Glauber model is un- 
derstood as the special case 5 = 0, F = 1.) 

The kink — > 3 kink processes are introduced via the 
exchange rate 

Wea;(Sj, Sj+l) = ^(1 - SjS.,+i) (2) 

This model (called NEKIM) , for negative values of S in 
eq.(l) shows a line of PC-transitions in the plane of the 
parameters {pex,S) 0. In NEKIMA H| the authors have 
extended the above model by introducing local symmetry 
breaking in the spin-flip rates of the -I- and — spins as fol- 
lows. Concerning the annihilation rates the prescription 
in is followed: 

wi+; ) = 1, w(-;++) = 0, (3) 

while further spin symmetry breaking is introduced in the 
diffusion part of the Glauber transition rate as follows. 
In calculating the transition rates 

p = w{-;+-)=w{-,-+) = T/2il-~6) (4) 

the Glauber form, eq.(Q, is used unchanged, while 
w{+; H — ) and w{+; —+) are allowed, to take smaller val- 
ues : 

p+ = w{+-+-) = w{+--+)<p. (5) 

In this way , by locally favoring the -|- spins, the effect of 
the other dynamically induced field arising from the pre- 
scription (cql^l) is counterbalanced. The spin-exchange 
part of the NEKIM model remains unchanged, eq.||2I). 
It is worth mentioning that the lifting of the strong re- 
striction in eq. © together with applying spin anisotropy 
p+ < r/2(l — S) has the same effect as a global magnetic 
field favoring -I- spins. 

For — 0, the absorbing states in the extreme situa- 
tion p+ = when diffusion-like spin- flipping maximally 
favors + spins, arc states with single frozen — spins like 

H h + H hH h++. By increasing p+ from zero, a 

slow random walk of these lonely — spins starts and by 
annihilating random walk only one of them survives and 
performs RW (see Fig^ . The all -I- and all — states are, 
of course, also absorbing. 




X 

FIG. 1: Space-time development of -I-' (white) and '-' (black) 
spins evolving from a random initial state for p+ = .1, pex ~ 0. 
Throughout the whole paper t is measured in units of Monte- 
Carlo sweeps. 



In ^2 the authors have studied the following global 
quantities for different values of p+ < p : the density 
of kinks as a function of time, starting from a random 
initial distribution of spins for pex — 

(6) 

and its asymptotic values for finite but small values of 

Pex 

PooiPex) Pex^ ■ (7) 

The results obtained within error of simulations, a = .5 
and f3 = 1.0, pointed to the presence of a N-BARW2 
transition. Finite size scaling behaviour was also exam- 
ined to find the other two bulk exponents, those of ^ , 
the correlation length and r, the characteristic time: 

i-Pex\r^e (8) 

where Z is the dynamical critical exponent. The expec- 
tation for a N-BARW2 transition at zero branching rate 
was justified by the values : i^j^ = 1.0 and Z = 2.0, which 
were found within error of simulations. We also found the 
expected phase diagram of a line of DP transitions in the 
{p+,Pex) plane (instead of the PC-line of NEKIM). 

III. CLUSTER BEHAVIOR AT AND BELOW 
THE N-BARW2 TRANSITION 

Spreading from a localized source at criticality is usu- 
ally described by the following three quantities 

p(^)-^-^n(^)-^^i?2(^)^^^ (9) 
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FIG. 2: Phase diagram of the NEKIMA model for 5 - 
r = .5. The absorbing phase is fuUy '— '. 
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FIG. 3; Cluster exponents for '-' cluster at p+ = .3, Pex ~ 

number of averages: 10^ "+" 



Number of MC steps: 10 
correspond, to 77, "x"-s to z/2, 



1 to S 



where n{t) denotes the mean population size, B?{t) is the 
mean square spreading of particles (here spins) about the 
origin and P{t) is the survival probability. In most cases 
these quantities are defined for particles, in the present 
case, however, like for studying compact directed perco- 
lation of an Ising chain ,15.] they will be used for spins. 

In the active phase the survival probability defines a 
further useful critical exponent (3 , [vn = Zvy^ ) 



P - r 



as 



(10) 



(11) 



In the present case the following parameter values of 
NEKIMA were used in the simulations: T = .5, 5 — 
— .565 {p = .39125). The phase diagram in the (A = 
, Pex) plane is shown on Fig. 2. The origin (0, 0) will 
be called "MDG-point" as at this point p+ = p and the 
model is the same as treated in [l3| (though the values 
of S, and F are different). The line Pex = is a line 
of compactness, as will be discussed in the next section. 
Also other details of the phase diagram will be explained 
later. 



5.10'^ MC steps and for averages over 10'' sam- 



pies. The local slopes 



<5(t) = 



\n[P{t)/P{t/m)] 
Inm 



(12) 



(and similarly for ri{t) and z{t)) as a function of 1/t are 
plotted, as usual in case of simulations for critically be- 
having quantities. In H12|l m > 1 is an arbitrary factor 
which we took to be equal to 5. The results obtained in 
different regions of the phase diagram, Fig|21 are sum- 
marised on Table I and Table II. As it is apparent from 
Table I, the cluster does not change its exponents by 
crossing the Pex = line. The '-' cluster's exponents, 
however, change abruptly. 

For the case p+ = .3 FigsI2| and 0] shows the local ex- 
ponent values, feall2|). for p^x — i.e. at the N-BARW2 
transition point and for = -02, i.e. in the active phase. 
Here and in most cases of our simulations the number of 
MC steps has been 5 x 10'^ with averaging over 2 x 10^ 
different runs. In some cases , however, much longer runs 
have also been carried out up to 10^ MC steps to corrob- 
orate these results, see Fig|5l 



In NEKIMA '+' and '-' spins are not symmetric, there- 
fore we have investigated two kinds of clusters. Namely, 
the development of the '-' cluster-seed was started from 
a wholly '+' environment at t=:0: ++++++++ - 
-|--|--t-+++++, while the cluster from a sea of '-' 
spins: 1 . We will call them '-' clus- 
ter and cluster, respectively. The simulations have 
been performed with several values of p+ and pex', for 



For comparison let us recall the well-known values for 
the above exponents in case of the compact directed per- 
colation point of the Domany-Kinzel cellular automaton 
[T^ . Dickman and Tretyakov |0 have given the results 
in this context as follows: 77 = 0, i5 = 1/2 and z = 1. 
( The same as for the Glauber-Ising model at 5 = 0, 
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TABLE I: Cluster critical exponents at and near the N-BARW2 transition point. The hyperscaling law rj + S = z/2 (see Section 
V.) is satisfied. 



exponents 


Pe. = 0, ' + ' 


Pe^ = '-' 


Pex / ' + ' 


Pex ^ '-' 


V 


1. 


0.0 


1.0 


1.0 


S 


0.0 


0.0 


0.0 


0.0 


z 


2.0 


0.0 


2.0 


2.0 



TABLE 11: Cluster critical exponents at and near the DP-line. For abbreviations see Fig. 2. The hyperscaling-law, valid for DP 
transitions, rj + 25 = z/2 (see Section V) is satisfied 



exponents 


on DP-line '+' 


on DP-fine '-' 


ABSO-phase '+' 


ABSO-phase '-' 


ACTIVE- phase '-f' 


ACTIVE-phase '-" 


V 


.31 


1.0 


exponential 


1.0 


1.0 


1.0 


5 


.16 


.0 


exponential 


0.0 


0.0 


0.0 


z 


1.26 


2.0 


exponential 


2.0 


2.0 


2.0 




0.003 



1/t 1/t 



FIG. 4: Cluster exponents for "-" cluster at p+ = .3, Pex = 
Number of MC steps: 10^, number of averages: 10^. "-|-"-es 
correspond to r\, "x"-s to z/2, *-s to 5 



FIG. 5: Cluster spreading exponents for cluster at param- 
eter values p+ = .3,pea; = .02 ( (5 = —.565, F = .5.) Number of 
MC steps: 10^, number of averages: 10^. Upper curve:r;e//, 
lower curve Ze///2. 



IV. BREAKING OF A SCALING LAW 



r = 1.0). It is of some interest to present the mea- 
sured cluster exponents at the origin of the phase dia- 
gram, Fig. 2, which is the (equivalent of the)MDG-point. 
Here we found for the cluster: ?/ = 0, (5 = 1/2, z ~\ 
while for the '-' cluster: r] = 1/2, i5 = 0, 2: = 1 (with the 
same accuracy as most of our results here i^max = 5.10^ 
MC steps)). These data are summarized on TABLE III. 
Because of the relatively low upper limit in time of most 
of our simulations as given above, the possibility of the 
presence of a log(i) correction at the MDG point can not 
be excluded. 



According to the previous section the result for the 
critical exponent of the mean square distance of spread- 
ing from the origin, z, is equal to 2.0 within error of 
numerical simulations. For the dynamical critical expo- 
nent the value Z = 2.0 was obtained, in the whole regime 
(P4. values ) of the N-BARW2 transition. 

On the other hand, at BARW-type transitions, like DP 
and PC transitions, the following scaling law connects the 
above two critical exponents: 

z = 2/Z (13) 

This relation is usually quoted as a consequence of dy- 
namical scaling. Using the above cited results, however. 
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TABLE III: Cluster critical exponents in case of the Glauber-Ising and MDG parameter values. The hyperscaling law rj+S = z/2 
is satisfied. 



exponents 


Glauber-Ising , '+' 


Glauber-Ising, '-' 


MDG, '+' 


MDG '-' 


V 








0.0 


0.5 


S 


1/2 


1/2 


0.5 


0.0 


z 


1 


1 


1.0 


1.0 



ea. (|13|l is broken. The possibility of breaking this scal- 
ing law is actually due to the circumstance, that the 
N-BARW2 transition point lies at the zero value of the 
branching probability, Pea; — and there is no absorbing 
phase, with exponentially decreasing space- and time de- 
pendences. To support this point let us recall the way 
Mendes et al. ^3 derived the relation H13|l . 

They started from the general expression for the den- 
sity of particles (kinks) at space-point r in the absorbing 
phase A < ( here A denotes the deviation from the 
critical point)and at large fixed value of t (for d = 1): 



(14) 



In the absorbing phase the function p{r, t) is expected 
to decrease exponentially as p{r,t) ~ exp— r/^, where 
^ ~ A""^^ . This form implies for F{u, v) (with v < 0) 
the form 



(15) 



where C > is constant. For to be time-independent 
the scaling law is required: 



2 
Z 



(16) 



This scaling law which is not fulfilled in the presently 
discussed model. Moreover, the bulk quantity, the time 
dependent kink density 

Pit) ^ (17) 

and the expression obtainable from ea. H14|) 

p{t) - (18) 

are also in conflict. Namely, while all the simulations 
have resulted in a = .5 within error and this is in agree- 
ment with the scaling law a = according to the val- 
ues given in Table I the exponent in eg. 1)18(1 is zero, again 
within the error of simulations. (It is to be noted, that 
this conflict is no more present concerning the exponent 
values at p^^ ^ 0, where a = and 77 — z/2 = 0, as well.) 
The apparent contradiction, however, is resolved again 
by recalling that cluster exponents and bulk-exponents 
are allowed to be not connected by a scaling law. 



V. HYPERSCALING 

The generalized hyperscaling law (l^j was developed 
for systems with multiple absorbing configurations and 
reads: 



2 1 



/3' 



S + 2t] ^ dz 



(19) 



where (3 is defined for the active phase, eq.(ll). 

The derivation of (|19|) goes along the following lines . 
It starts with ea. l|14|) for p{r,t) and with the expressions 



P{t) r^t-^'S>{At^/''«) 



A^' ,l3'=5v 



(20) 



(21) 



for the survival probability. Since the stationary distri- 
bution is unique 



p{x,t) P^A^ ^ l^P+P 



(22) 



as t ^ 00. Hence -^(0, y) ~ yl^+l^ which entails ea. (|19|l . 

In case of the DP transition (along the DP line of Fig. 2) 
(3 = /3 as it is well known, and thus ea. ((T^ gives 



25 + 7? = z/2 



(23) 



For the N-BARW2 transition, however, ea. H19|l does 
not apply as ea. lfT^ . according to the previous section, 
is not an appropriate starting point. 

To deduce the hyperscaling law valid for this case there 
are several possible ways of arguing. It is possible to 
enlarge the parameter space of our model: we can think 
of a third direction in the parameter space, approaching 
from where the transition turns out to be of first order. 
For this aim one can introduce a 'magnetic field' into 
the system by changing the annihilation probability as 

^(-f; ) = \ — h. In this direction j3h = and thus 

ea. ((T^ gives (d=l): 



(24) 



This law is satisfied for all the clusters investigated, in- 
cluding those at the MDG point. Even for pex ^ we can 
still think of each point as being a first order transition 
point with /3 = in the /i-direction and the same consid- 
erations apply as above. Thus on the basis of the results 
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presented now, the conclusion to be drawn is that hy- 
perscahng is generically satisfied in the whole N-BARW2 
phase of the NEKIMA model. 

Looking at the problem from a different point of view, 
however, it is really not necessary to introduce the above 
auxiliary magnetic field. Namely, one can simply make 
the observation, that all the clusters investigated on the 
Pex = line are compact and from this fact Eq|23Jl fol- 
lows for the hyperscaling law [l5| . 

Ea. 1)24(1 is known as the hyperscaling law for compact 
clusters. By definition S+rj is the exponent which charac- 
terizes the average population in surviving trials and the 
radius of such a cluster grows as Rt ~ t^^^. 5 -I- 77 = c?| is 
simply the scaling law for the volume of a d-dimensional 
sphere of radius Rt ^3 . 

As a matter of fact, the '-' clusters are compact with 
5 — also at p^x 7^ 0, and even in the DP-region of Fig. 2. 
This is not true, however, for the '+' cluster in the DP- 
region, which follows normal DP-cluster behavior (see 
TABLE II) with the corresponding DP hyperscaling law 
eq. (|23|l (It is worth noting that whenever (5 = 0, and this 
fails only for the cluster in the DP-region, the CDP 
and DP hyperscaling laws do not differ). 



VI. REENTRANT PHASE DIAGRAM, 
CLUSTER MF CALCULATIONS 



In the original NEKIM model at 5 > no transition 
occurs, while for negative values of this parameter PC 
transition takes place. The spin asymmetry of NEKIMA 
changes the character of the transition into DP and this 
appears also for (5 > Here we have chosen for our simu- 
lations and for the cluster mean-field approximation cal- 
culations a fixed negative value of 5. Our aim has been 
to explore some possible reminiscence of the PC transi- 
tion. At the chosen parameter values 5 = — .565,r — .5 
in NEKIM the PC transition occurs at p^x — -12 Turn- 
ing to NEKIMA, at the same values of (5, F our simu- 
lations show that the transition point (which is DP, of 
course) shifts to Pex — -51. The absorbing phase below 
this point is all For letting p^ < p the DP-line starts 
tangentially upon increasing p^x from and exhibits a 
reentrant property. It ends up at Pex = .51 tangentially. 
The regression takes place at Pex ~ -12, see Fig. 2, most 
probably a remnant of the transition point of the corre- 
sponding PC transition. This turning point, however, is 
also of DP character as can be expected. 

Dynamical cluster mean-field approximations have 
been introduced for nonequilibrium models by l20l| . 
The master equations for = 1 — 7 block probabilities 
were set up 



dPN{{si\) 
dt 



f{PN{{s^})) 



(25) 
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0.08 



0.06 - 



0.04 



0.02 




ACTIVE 



ABSORBING 



where site variables may take values: Si — ±1. Taking 
into account spatial reflection symmetries of PN{{si}) 
this involves 72 independent variables in case oi N = 



0.4 0.395 0.39 0.385 0.38 0.375 0.37 
P. 

FIG. 6: Steady state density in A = 6 (bullets) and A = 7 
(boxes) level approximation. Lines connecting symbols are 
shown for guidance of eye only. 



7. The equations were solved numerically for the 
^^"ot"'"^"* ~ ^ steady state condition, for different Pex 
and p+ values and the Pkipo) kink density was expressed 
by P/v({sj}). The reentrant behavior could not be ob- 
served for N < Q clusters. The results for A^ = 6, 7 are 
shown on Fig. |S| A slow shift towards lower p_|_ values, 
which agrees with the simulations can be observed. 



VII. DISCUSSION 

We have investigated a one-dimensional nonequilib- 
rium spin model (NEKIMA) exhibiting strong spin- 
asymmetry. In the plane of two of the parameters of 
NEKIMA (the kink-branching parameter and a spin- 
asymmetry parameter) the phase diagram is as follows: 
besides a reentrant DP line the NBARW-2 transition oc- 
curs at zero branching rate. Due to the asymmetries, 
and '-" spin clusters behave differently. By investigating 
their development we conclude that generic power-law 
behavior characterizes the cluster behavior at and in the 
vicinity of the NBARW-2 transition. 

The critical cluster exponents obtained satisfy the con- 
straints on critical exponents in general: 1. (5 > , 2. 
1 < z < 2. The critical exponent 5 has been found to 
be zero. The hyperscaling law is satisfied in the form 
known for compact directed percolation, and indeed, the 
N-BARW2 clusters are compact. 

In a different problem Cafiero et al. |2l| have re- 
ported cluster exponents similar to the ones found here. 
These authors studied how disorder affects the critical 
behaviour of DP-like systems. Already in the eighties 
Noest showed that quenched disorder changes their 
behaviour in c? < 4 and demonstrated that when d = 1 
a generic scale invariance can be observed. In 0| it 
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was shown that deep in the active phase rj = 1, 5 = 
and z = 2 for the model they considered. As we have 
found also generic scale invariance and the same expo- 
nents, in the active phase of our model and even in the 
region which is the active phase of the DP line of our 
phase diagram, the question arises whether the similar- 
ity is fortuitous or not. Whether the slowly diffusing 
'-' clusters of NEKIMA distributed randomly in the x- 
direction can play a role similar to quenched impurities 
e.g. in the original NEKIM model is a question for future 
investigations. 
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